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ABSTRACT 


Extraterrestrial  radio  waves  from  sources  of  small  angular  size  can  be  used 
to  synchronize  clocks  at  two  receiving  sites  to  a  high  degree  of  precision. 
Radio  star  signals,  together  with  time  identification,  are  recorded  on  mag¬ 
netic  tape.  Tapes  are  subsequently  processed  on  a  digital  computer  to  obtain 
the  synchronization  error  between  the  two  clocks  at  the  time  of  observation. 
The  ultimate  limit  to  the  precision  of  the  method  is  the  knowledge  of  the  rel¬ 
ative  signal  delay  in  the  atmosphere  and  ionosphere. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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VERY  LONG  BASELINE  INTERFEROMETRY 
AS  A  MEANS  OF  WORLDWIDE  TIME  SYNCHRONIZATION 


I.  INTRODUCTION 

1-3 

Interferometers  with  baselines  up  to  7,000km  have  been  used  to  measure  the  angular 
sizes  of  extraterrestrial  radio  sources.  These  interferometers  use  independent  frequency 
standards  of  a  high  spectral  purity  to  convert  radio  star  signals  to  video.  The  video  signals 
from  each  element  of  the  interferometer  are  encoded,  together  with  time  identification,  and  re¬ 
corded  on  magnetic  tape.  The  tapes  are  subsequently  processed  on  a  digital  computer  to  obtain 
information  about  the  angular  size  of  the  radio  source. 

H.  PRINCIPLES  OF  OPERATION 

If  the  angular  size  of  a  radio  source  at  wavelength  X  is  smaller  than  x/l)  radians,  the  elec¬ 
tric  field  time  functions  sampled  at  a  distance  D  apart  along  the  wavefront  are  identical.  Thus, 
the  voltage  time  functions  received  by  two  antennas  from  a  small-diameter  source  are  just 


Fig.  1.  Long  baseline  interferometer. 


delayed  replicas  of  each  other.  This  is  illustrated  in  Fig.  1.  The  delay  r.^  is  determined  by 
the  geometry  shown  in  Fig.  2.  Neglecting  delays  in  the  atmosphere  and  ionosphere,  r.^  is 
simply 

D/C  cos  9 

where  D  is  the  length  of  the  baseline,  C  is  the  velocity  of  propagation,  and  G  is  the  angle  be¬ 
tween  the  line  to  the  source  (assumed  at  infinity)  and  the  line  between  sites.  In  celestial 
coordinates, 

Tint  =  T'  s*n  +  cos  cos  COS  <LS  ~~  )] 

where  Lg  and  are  the  hour  angle  and  declination  of  the  source,  and  Lg  and  6^  are  the  hour 
angle  and  declination  of  the  baseline. 
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Fig.  2.  Geometry  of  long  boseline  interferometer. 

The  geometric  delay  could  be  measured  by  cross  correlating  the  time  functions  known  to 
occur  at  specified  times.  Alternatively,  the  difference  in  time  between  clocks  used  to  identify 
time  function  samples  can  be  obtained  if  the  geometric  delay  is  already  known.  In  practice, 
the  time  functions  cannot  be  sampled  directly  but  must  be  converted  to  video  with  a  local  oscil¬ 
lator  and  mixer  chain. 

For  a  continuum  source  (a  source  whose  spectrum  is  nearly  flat  or  nonspectral),  the 

Fourier  transform  S  (cj)  of  the  cross-correlation  function  R  (r)  of  the  video  signals  (in  an 

xy  xy 

upper  sideband  system)  is 


Sxy(w)  =  A  exp  {i[<f>  +  Awt  +  (w  +  wQ)  (rint  +  )]} 

where 

=  local  oscillator  frequency  (radians/sec), 
cj  =  local  frequency, 

Acj  =  difference  in  local  oscillator  frequencies, 

r  =  time  error, 
e 

(p  =  instrumental  phase, 

cj  dr.  Vdt  =  fringe  rotation  rate, 
o  int  & 
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This  function  has  a  time- varying  phase  given  by  the  sum  of  the  local  oscillator  frequency  dif¬ 
ference  and  the  "fringe  rate."  This  can  be  compensated  by  suitable  local  oscillator  offsets  (to 
cancel  most  of  the  fringe  rate)  and  by  sine-cosine  multiplication  in  the  data  reduction  program. 
However,  a  fundamental  "coherent"  integration  time  T  is  imposed  by  requiring 

l<e^(t)>T|«l  , 

or  the  "phase  noise"  between  standards  must  be  small  enough  to  allow  coherent  integration. 
Figure  3  shows  how  long  a  coherent  integration  can  be  performed  at  1666  MHz  by  requiring  A <p 


Fig.  3.  Fractional  frequency  stability  and  phase 
noise  of  various  frequency  standards.  Data  were 
obtained  by  R.  Vessat,  Hewlett-Packard  Fre¬ 
quency  and  Time  Division,  November  1967. 


to  be  less  than  100  degrees.  The  cross-correlation  function  R  (t),  after  removal  of  fringe 
rotation  rate,  is 


R  ( r )  =  2A  cos  (p 
xy 


sinamax(T  +  Tint  +  Te> 

*(T  +  Tint  +  Te) 


—  2A  sin  (p 


[l-cosa,max(T  +  Tintf  Te)] 

*{T  +  Tint  +  Te> 


where  cj  is  the  bandwidth  (radians/sec)  of  the  video  recorded.  The  above  shows  that  the 
max  ' 

cross  correlation  is  zero  except  around  r  =  —  (r.^  +  re).  Figure  4  shows  sine  and  cosine  com¬ 
ponents  of  the  cross-correlation  function. 

To  reduce  the  information  rate,  and  for  ease  of  computation,  the  video  signals  are  first 
bandlimited  and  then  infinitely  clipped  and  sampled  at  a  rate  equal  to  twice  the  bandwidth.  Thus, 
the  video  data  are  in  binary  form  (video  ^  0  corresponding  to  logical  1,  video  <  0  corresponding 
to  logical  0)  and  are  easily  digitally  recorded*  in  the  form  of  bit  strings.  The  effect  of  infinite 

4 

clipping  is  to  preserve  zero-crossing  information,  and  a  theorem  by  Van  Vleck, 


*At  present,  time  identification  is  achieved  by  recording  the  first  bit  af  each  data  block  at  a  precisely  known 
time  relative  ta  the  station  clack. 
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SIN  CORRELATION  COS  CORRELATION 


Fig.  4.  Cross-correlation  function  components  observed  on  radio  source  3C273 
with  receiver  bandwidth  of  120  Hz. 
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Pxy(r)  =  sin^/Z  pc 


(r )] 


xy 


allows  the  normalized  cross-correlation  function  of  the  video  signals  p  (t)  to  be  estimated  from 

xy 

the  cross  correlation  of  the  bit  strings  pc  (r).  Figure  5  shows  a  block  diagram  of  a  long  base- 

xy 

line  interferometer.  Figure  6  shows  the  amplitude  of  the  S  (w)  for  different  trial  fringe  rates. 

xy 


Fig.  5.  Long  baseline  interferometer  terminal. 


Fig.  6.  "Fringe  rate"  spectrum  of  3C273  observed  with  coherent  integration  time 
of  160  seconds  compared  with  theoretical  |(sin  x/x)|  function  which  would  be  ob¬ 
served  in  the  absence  of  phase  noise. 
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WIDTH  OF  MAINLOBE  80 nsec 


Fig.  7.  Theoretical  delay  resolution  function  for  multiwindow  system. 


Fig.  8.  Sidelobe  structure  of  function  shown  in  Fig.  7. 
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HI.  DELAY  RESOLUTION 


From  the  expression  for  the  cross-correlation  function,  it  is  clear  that  the  delay  resolu¬ 
tion  is  of  the  order  of  (27r/wmax)>  since  this  is  approximately  the  half- power  width  of  the  mag¬ 
nitude  of  the  cross-correlation  function.  If  the  signal- to- noise  ratio  is  good  enough,  resolution 
should  be  up  to  one  order  of  magnitude  better.  However,  if  recorder  bandwidth  is  limited,  the 
delay  resolution  can  be  improved  by  sampling  and  multiplexing  several  narrowband  video  win¬ 
dows  within  a  wide  band.  When  this  procedure  is  adopted,  the  delay  resolution  function  (cross¬ 
correlation  functions  convolved  with  expected  functions  and  summed)  is  shown  in  Figs.  7  and  8 
for  one  proposed  multiwindow  scheme. 


IV.  SYSTEMATIC  ERRORS 

The  long  baseline  interferometer  as  a  precise  method  of  time  synchronization  is  limited  by 
the  precision  with  which  the  geometric  delay  can  be  computed.  The  computation  of  geometric 
delay  involves  a  knowledge  of  source  positions,  baseline  vector,  and  sidereal  time.  Without 
first  using  the  interferometer  as  a  survey  instrument,  the  geometric  delay  can  only  be  computed 
reliably  to  within  2  0  nanoseconds.  For  baselines  within  a  continent,  it  might  be  possible  to  re¬ 
duce  errors  below  2  0  nanoseconds.  With  errors  this  large,  it  is  not  important  to  take  careful 
account  of  ionospheric  and  atmospheric  delays  which  should  not  amount  to  more  than  2  0  nano¬ 
seconds  at  1.6  GHz  except  at  low  elevations.  However,  a  careful  use  of  a  series  of  observations 
aimed  at  determining  both  time  synchronization  and  baseline  vector  should  result  in  an  order  of 
magnitude  improvement  if  the  ionospheric  delay  can  be  correctly  predicted. 


V.  SIGNAL-TONOISE  RATIOS 

Up  to  this  point,  we  have  considered  the  video  signals  to  be  due  only  to  the  radio  source. 
With  moderate  antenna  sizes  and  conventional  radiometric  receivers,  the  noise  far  exceeds  the 
source  signal.  However,  the  noise  signals  are  uncorrelated  and,  provided  that  enough  bits  are 
processed,  it  is  possible  to  observe  very  weak  signals.  The  signal- to- noise  ratio  (S/N)  after 
coherent  integration  is 


S_ 

N 


T  T 
a  a 


T 


2  \Tbt 

y 


where  the  first  factor  is  the  geometric  mean  of  the  signal- to- noise  ratio  in  the  two  systems 
before  data  processing,  and  the  second  factor  is  l/s/T  times  the  square  root  of  number  of  bits 
processed.  Table  I  gives  the  signal-to- noise  ratios  for  various  antenna  combinations  for  100 

7 

seconds  integration  or  7.2  x  10  bits  processed.  If  a  signal- to- noise  ratio  of  4  is  regarded  as 
adequate,  a  combination  on  the  NRAO  140-foot  and  any  2  0- foot  antenna  with  100°K  system  tem¬ 
perature  could  yield  this  desired  signal- to-noise  ratio  with  only  about  10^  bits  processed. 

106  bits  could  be  transmitted  over  a  normal  HF  radio  channel  in  about  5  minutes. 


7 


TABLE  I 

SIGNAL-TO-NOISE  RATIO  AT  1.6  GHz  ON  3C273 
(100  Seconds  Integration  or  7  X  10^  Bits  Processed) 


Antenno  1 

Antenno  2 

"A,  i 

(1VVAnt  2 

S/N* 

Hoys  lock  120-foot 

NRAO  140-foot 

1/100 

1/23 

100 

Hoystock  120-foot 

Any  60-foot  ontenno 

1/100 

1/200 

40 

NRAO  140-foot 

Any  20-foot  ontenno 
with  100°K  system 

1/25 

1 

(49  X  25) 

40 

*S/N  bosed  on  14  flux-units  unresolved  from  3C273  ot  1.6  GHz. 


Other  Strong  Sources  in  Unresolved  Flux-Units 


CTA21 

6 

3C287 

3 

3C84 

3 

3C286 

3 

NRAO  140 

3 

3C309. 1 

3 

NRAO  150 

4 

1510-08 

3 

3C19 

3 

3C345 

5 

3C120 

3 

NRAO  530 

4 

3C138 

3 

3C380 

3 

3047 

4 

3C418 

4 

1127-14 

6 

2127+04 

4 

1148-00 

3 

2203-18 

3 

3C268. 1 

3 

3C446 

3 

3C273 

14 

CTA102 

6 

3C279 

8 

3C454. 3 
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VI.  CONCLUSION 

The  very  long  baseline  interferometry  technique  should  be  able  to  provide  worldwide  time 
synchronization  to  0.05 psec  for  time  standards  located  at  the  world's  major  radio  astronomical 
observatories  and  surveillance  radar  establishments.  At  first,  data  could  be  recorded  on 
standard  computer  tape  and  processed  on  general-purpose  computers.  However,  if  synchroniza¬ 
tion  were  required  on  a  close  to  real-time  basis,  HF  radio  links  could  provide  data  links  to  a 
computer. 
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